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We introduce state-independent, non-perturbative Hamiltonian quantum speed limits for population leakage 
and fidelity loss, for a gapped open system interacting with a reservoir. These results hold in the presence of 
initial correlations between the system and the reservoir, under the sole assumption that their interaction and 
its commutator with the reservoir Hamiltonian are norm-hounded. The reservoir need not he thermal and can 
be time-dependent. We study the significance of energy mismatch between the system and the local degrees of 
freedom of the reservoir which directly interact with the system. We demonstrate that, in general, hy increasing 
the system gap we may reduce this energy mismatch, and consequently drive the system and the reservoir into 
resonance, which can accelerate fidelity loss, irrespective of the thermal properties or state of the reservoir. This 
implies that quantum error suppression strategies based on increasing the gap are not uniformly beneficial. Our 
speed limits also yield an elementary lower bound on the relaxation time of spin systems. 


Quantum speed limits (QSLs) answer the fundamental 
question of how fast a quantum system can evolve, and have 
numerous applications, e.g., in quantum computation, control, 
and metrology. Traditionally, they characterize the minimum 
amount of time required for a quantum state of a closed quan¬ 
tum system to evolve to an orthogonal state. Mandelstam & 
Tamm (MT) Ul first showed that this time is lower bounded 
by the inverse of the standard deviation of the Hamiltonian. 
Margolus & Levitin (ML) ||2l gave a different bound involv¬ 
ing the inverse of the mean of the Hamiltonian, and the two 
bounds were subsequently unified |[3. These results led to 
numerous applications and extensions which go beyond the 
traditional QSLs, and consider, e.g., the minimum time for 
optimal control, or for implementing a unitary gate in quan¬ 
tum computation Il4l 420l . 

In this Letter we are concerned with QSLs for open quan¬ 
tum systems ED, a question that has attracted significant re¬ 
cent attention Il22l4^ . While earlier work focused on gener¬ 
alizing the MT or ML-bounds to the open system setting, here 
we present state-independent, non-perturbative Hamiltonian 
QSL bounds for population leakage and fidelity loss, for a 
gapped open system interacting with a reservoir. The assump¬ 
tions behind the results we present here are also different and 
independent from those behind previous such bounds. First, 
we make the (often natural) assumption that the system’s ini¬ 
tial state is restricted to an energy sector which is separated 
from the rest of the spectrum by a nonzero gap AE, e.g., the 
ground subspace in various quantum information processing 
applications. Second, our bounds are independent of the state 
of the system or reservoir, and in particular, remain valid in 
the presence of initial correlations between the system and 
the reservoir. Third, our bounds are obtained purely at the 
Hamiltonian level. Thus, unlike most other open system QSL 
bounds II 22 I - E 5 I . we do not assume that the system’s evolu¬ 
tion is governed by a master equation or a completely positive 
channel. 

Our key result is given in Eqs. @ and Q below, and com¬ 
prises fundamental QSL bounds on decoherence and leakage 
times, expressed in terms of AE and the bounded norms of the 


interaction Hamiltonian and its commutator with the reservoir 
Hamiltonian, without assuming that the reservoir Hamiltonian 
is norm-bounded. Note that applying the traditional closed 
system QSL bounds to the system and reservoir together in 
general yields bounds which are rather loose and independent 
of the gap AE E^ . 

Given the very general assumptions behind our QSLs, they 
have a wide range of applicability similar to the previously 
known QSLs, including relaxation in many-body spin systems 
and limitations of control via a remote controller. The primary 
application on which we focus is quantum error suppression, 
where the goal is to slow down the loss of fidelity relative 
to some desired system state, e.g., in the context of quantum 
information processing tasks isiiiii. A common strategy 
to achieve fidelity enhancements is to use or generate energy 
gaps (e.g., II 29 H 35 ]). Therefore, after deriving our QSLs we 
study the dependence of the speed of decoherence and leakage 
on AE. This enables us to find a general lower bound on the 
timescale for leakage. As expected, we find that in the AE —t 
00 limit the probability of leakage at any finite time goes to 
zero, and moreover, that if the error detection condition ll^ 
holds then in this limit the state retains its fidelity and remains 
unaffected by the environment. However, we demonstrate that 
such protection is not guaranteed when AE is finite. Namely, 
by analyzing a spin system model, we show that increasing 
the gap can in fact accelerate fidelity loss and decoherence, 
essentially because of a resonance between the system and 
the reservoir. This means that protection via increasing energy 
gaps can be counterproductive 137). 

Technical results .—Consider a system S coupled to a reser¬ 
voir R with the total Hamiltonian iTtot(i) = Hs+ H^{t) + Hi 
where [Hs, HR{t)] = 0 and the interaction satisfies ||iTi|| < 
00 (we use the operator norm || • |j, i.e., the largest singu¬ 
lar value; we also use h = 1 units throughout). An impor¬ 
tant class of examples are spin-bath models We denote 
the time-dependent joint system-reservoir state evolving un¬ 
der iTtot(f) by psR(f) the reduced state of the system at 
time t by p{t) — TrR[psR(i)]- Let C be the subspace of the 
system Hilbert space spanned by the eigenstates of Hs whose 


2 


energies lie in the interval I C K, which includes at least 
one eigenvalue of Hs. Let Pc be the projector onto C, and 
Qc = I — Pc he the projector onto the orthogonal subspace 
C-^. Thus [Pc,Hs] = 0. Let SE denote the energy spread 
in C, i.e., the difference between the minimum and maximum 
eigenvalues of Hs in I. Let AE denote the gap between the 
energy levels of Hs inside and outside C (i.e., if Ai and A 2 are 
two distinct eigenvalues of Hs such that Xi G X but A 2 ^ X, 
then |Ai — A 2 I > AE). Initially we assume AE > 2||i7i||, 
which guarantees that there is a separation between the sys¬ 
tem energies inside X and the rest of the spectrum, even in the 
presence of the interaction. This simplification is relevant be¬ 
cause we are mostly interested in the large AE limit. Later, 
when we arrive at Eq. o, we present the general form of 
the result which relaxes this condition, and results in tighter 
bounds for small t, even when AE < 2||ffi||. Before we in¬ 
troduce our bounds, we define an important inverse timescale 
for open system dynamics, that will make repeated appear¬ 
ances: 


n{t) = 


2||[gi,gR(t)]|| 

AE-2\\Hi\\ 


( 1 ) 


We proceed to present and interpret our main results. All our 
results are given rigorous proofs in the Supplementary Mate¬ 
rial (SM) ll^ . Unless stated otherwise, throughout we as¬ 
sume that the system state is initialized in C, i.e., p(0) = 

Pcp{h)Pc- 

Leakage .—Leakage is the process whereby the system state 
develops support in C'*', which we quantify in terms of the 
leakage probability pieak(i) = Tr[p(f)Qc]- Our first general 
result is an upper bound on pieak(f)j proved in the SM Il39l : 


Pleak(i) < 


(mi\ 


V AE 


ds U(s) 


AB- 


0 


( 2 ) 


To explain this bound, note that the terms ||Lfi||/Aiii and 
/p ds U(s) correspond to two different sources of leakage: 
||Lfi||/Ai? determines how much C is rotated by the interac¬ 
tion Hi. The rotated eigenstates of the perturbed Hamiltonian 
can cause leakage relative to the eigenstates of the original 
Hamiltonian. Of course, this also happens in the closed sys¬ 
tems, and this is why this term does not vanish for H^{t) = 0, 
where the total system Hamiltonian becomes Hs + Hi. Since 
||Lfi||/Ai? is time-independent, it remains small and insignifi¬ 
cant in the limit where the gap is large. The term ds r2(s) is 
more interesting. In particular, in the case of time-independent 
Hn{t) = Hu, where the total energy of the system and reser¬ 
voir is a conserved quantity, || [H^, iTi]|| can be interpreted as 
the maximum rate of change of energy of reservoir. Then, 
in the special case where C is the bottom (top) energy sec¬ 
tor, 0“^ can be interpreted as the minimum time the reservoir 
needs to transfer (absorb) the required energy to move the sys¬ 
tem from C to (see the SM 09l ). 

Fidelity .—We compare the instantaneous “actual” state 
pit) and the “ideal” system state pid(f) = 
using their Uhlmann fidelity HOI |4T] F[pit), puit)] = 
\\\/pit)\/Pidit)\\i (II ■ 111 is the trace norm) and their Bu¬ 


res angle 0(f) = arccos (F[p(f), pid(f)]), a generalization to 
mixed states of the angle between two pure states ll42l . Let 
Pq = Pc ^ Ir. We define the induced splitting by Hi on C as 


IS(Po^^i-Po) = ininjfRGHerm(WR)||-Po-ffi^o — Pc® .?^r|| > (3) 


where the minimization is over the Hermitian operators act¬ 
ing on the reservoir Hilbert space Hr. This quantity can be 
interpreted as the strength of the effective interaction between 
the code subspace and the reservoir in the lowest order of per¬ 
turbation theory. It exists because the reservoir can couple to 
different states in the subspace C in different ways and this 
generally leads to decoherence, or, in special cases, to a mod¬ 
ification of the system Hamiltonian, a potentially beneficial 
effect ll4^ (see the SM lf39l l. This term can be nonzero only 
when C is at least two-dimensional. We can now state our 
second general result, an infidelity upper-bound: 


sin^ = ^yi-F[p(t),p.d(t)] < M (4) 


-f f ds H(s) t 

Jo 


ISjPoHiPo) 

2 


2\\Hi\\iSE+\\Hi\\y 

AE 


Related bounds have been obtained in Ref. Il34l . While 
bound Q holds for AE > 2||iTi|| and states initialized in 
C, our third general result is a simple universal QSL bound 
which does not require either one of these assumptions: 


. 0(f) f(Amax-Amin) f||Lfl| 

sm —— < -;- < —-— 

9 - A - 9 


(5) 


where Amax and Amin are the maximum and minimum eigen¬ 
values of Hi, respectively. This bound formalizes the stan¬ 
dard intuition that the minimum relaxation time of an inter¬ 
acting system is determined by the inverse of the couplings. 
However, as we will show in an explicit example, our QSL 
bounds (|^ and Q can lead to much stronger bounds on the 
relaxation time. 


Quantum speed limits .—The bounds we have presented di¬ 
rectly lead to QSLs on open-system quantum evolution, as we 
show next. For simplicity, in the following we assume that 
HRit) = Hr. 

Let denote the smallest time at which the probability 
of leakage from C exceeds a constant threshold po G (0,1). 
Then, it follows from bound (|^ that in the large-gap limit 
(i.e., ||7Ti||/Ai? <C Pq ) this timescale is lower-bounded by 
2 ||[.^^r]|| Po^^' ^ different lower bound on 

using bound together with the fact that U[p(f), Pid(f)] < 
\/l — Pieak(i) (see the SM ll^ l. Let Tmin be the small¬ 
est time at which F[p(f), pid(f)] drops below the threshold 
(1 — po)^^^ for an arbitrary initial state. This threshold con¬ 
vention guarantees > r^in- Then bound (|^ implies 
Tmin > c(po)IILfi|i“\ where c(po) = [2{1 - (1-po)^/^}]^/^, 
and hence 

TiL > max|c(po)||iTi|r^ , 2\\[Hi^Hr]\\ } ' 
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Similarly, we can define to be the smallest time at which 
F[p{t), Pid(t)] drops below the threshold (1 —For this 
choice of threshold we always have t£j < If we further 
assume the same large-gap limit and also that IS(Po-ffiFb) = 
0 and 5E = 0, which is a relevant assumption in the context 
of error suppression, we find using bound Q that 


T-fid > c(po) max<^ \\Hi 


1-1 


AE 


4(li[iFi,iFR]|| 


\Hi\n 


(7) 


Equations (|^ and Q constitute our key new QSL bounds. It 
follows from the definitions of the various timescales we have 
introduced, together with our result in the bound (|^, that 


eL > Tfid^T-min > c(po)||i7l|| ^ • (8) 


The above bounds on and Tmin are all first-order in 

On the other hand, any master equation derived 
under the Born-Markov approximation (BMA) is necessarily 
second-order in the coupling strength ED. Therefore, these 
QSL time scales, or more generally any open-system behav¬ 
ior which occurs on a timescale of order such as the 

resonance phenomenon discussed below, cannot be described 
under the BMA. 

Quantum error suppression .—One of the main applications 
of these bounds is in the context of quantum error suppres¬ 
sion. C is then the code subspace and one is usually interested 
in the case where it is a degenerate eigensubspace of (i.e., 
5E = 0). In this case pid(f) = p(0), whence F[p(f), pid(f)] 
is simply the fidelity between the initial state and the state at 
time t. The fidelity can degrade even if the gap is large com¬ 
pared to the interaction, i.e., if ||Tfi|j/Ai? < e ^ 1. In this 
limit bound Q implies that the rate of fidelity loss is upper 
bounded by O'(f) = 2y20(f)-bv/2IS(Po-ffiPo)+0(e)||-ffi||. 
This result has a simple interpretation: fidelity loss can hap¬ 
pen either because of leakage, whose speed is bounded by 
0 (t), or because of the effect of the reservoir on C, whose 
strength is given by IS^PqHiPq). In the limit AE oo the 
rate of fidelity loss is determined just by the induced splitting 
lS{PoHiPo), and if this quantity vanishes then for any finite 
time t, F[plt),pid{t)] 1. 

Therefore, the special case where IS{PoHiPq) = 0 is par¬ 
ticularly important for error suppression. To illuminate it, 
consider the decomposition of the interaction term as Fli — 
® where {S'"} and {B°‘} are, respectively, inde¬ 
pendent system and reservoir operators. Then, using Eq. •ID, 
we find that lS{PoHiPo) = 0 iff PcS^^Pc oc Pc for all S". 
This is also known as the quantum error detection (QED) con¬ 
dition ll^[^ . Thus the induced splitting quantifies the devi¬ 
ation from the QED conditions. \S{PqF[\Pq) — 0 can be the 
result of symmetries of the interaction as in a decoherence- 
free subspace ll44l |45]| . or it can be engineered using QED 
codes (e.g., EH)- Using bound 0 we can go beyond this 
special case and study the effectiveness of a particular error 
suppressing scheme in the case where the perfect QED condi¬ 
tion does not hold (see also Ref. Il46l i. 

Role of the reservoir and system parameters .—One of the 
interesting aspects of bounds (|^-Q is that they are indepen¬ 


dent of the reservoir state, and the reservoir Hamiltonian en¬ 
ters only via || [iJi, iJR(f)] |j. This means that even if the reser¬ 
voir is infinitely large and ||iTR(f)|| or \\dHj^{f)/dt\\ are un¬ 
bounded, as long as i/R]|| remains small and bounded, 
the leakage can be a slow process, depending on the ra¬ 
tio Ai5/|| [Tfi, iTR]||. This happens, in particular, when the 
interaction with the reservoir is quasilocal, i.e., the system 
degrees of freedom (DOEs) interact only weakly with the 
distant reservoir DOEs. To be concrete, consider the de¬ 
compositions Ell = where each term acts 

non-trivially only on a local DOE i in the reservoir. Then 

muH^m < 2E.eR 11 ^ 1*^11 where is 

the sum of all the terms in which act non-trivially 

on DOE i. In many physical scenarios this sum, and hence 
II[Tfi,TfR(f)]||, is bounded and small while ||i/R(t)|| is un¬ 
bounded and contains long-range interactions. E.g., the reser¬ 
voir may contain bosonic DOEs, for which ||i/R(f)|| is infi¬ 
nite. But, as long as these bosonic DOEs do not directly inter¬ 
act with either the system or the DOEs which directly couple 
to the system (i.e., those with 0), || [i/j, 7TR(f)]|| can 

be small. This remains true even if information propagates 
arbitrarily fast through the reservoir and the Lieb-Robinson 
velocity m is unbounded. 

On the other hand, if || [Hi, Hn{t)] |j is large and comparable 
to A£'||iJi||, then our QSL bounds suggest that the timescales 
for fidelity loss and leakage error can be as small as 
even in the large AE limit. As we explicitly show below, the 
bounds are attainable when || [Hi, iTR(f)] |1 — Ai?||iTi||. 

It is also interesting to note that bounds (|^-Q are inde¬ 
pendent of the state of the reservoir. This implies that even in 
the limit of infinitely high temperature T, leakage can still be 
a very slow process, depending on the ratio AE/[\ [Hi, iTRUj. 
This is a consequence of the assumption that both ||i7i|| and 
II [Hi, TIr] II are bounded, and it does not hold, e.g., in the case 
of bosonic reservoirs with the standard spin-boson coupling to 
the system. On the other hand, even at T = 0 all information 
in the system state can be erased by the reservoir in a time 
of order ||iJi||”^, the shortest timescale over which the reser¬ 
voir can have any influence on the system. Thus, the time it 
takes the reservoir to affect the evolution of the system is not 
necessarily related to T. 

A model of resonance .—To study the dependence of the 
time scales for leakage and fidelity loss on the system and 
reservoir parameters, we present an illustrative example. This 
is a simple model that exhibits the phenomenon of resonance 
between the system and reservoir, and is relevant, e.g., also in 
the context of state transfer via spin chains ||48]| . The system 
is a single qubit {k = 1) with Hamiltonian AEiaf/2, and 
gap AEi. The reservoir can have an arbitrarily large number 
of DOEs and may contain bosonic modes. The only assump¬ 
tions we make about the structure of the reservoir are (i) the 
only reservoir DOE which directly couples to the system is an¬ 
other qubit (k = 2), and (ii) the interaction between the reser¬ 
voir qubit and other reservoir DOEs, denoted by /i 2 ,rest (t), is 
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bounded. The total Hamiltonian is 

^ AT/ 

^ +J*?! •a2+/l2,rest(i)+i^rest(0 ) (9) 


where iTrest(i) is an arbitrary Hamiltonian that acts trivially on 
qubits 1 and 2. 

The system qubit is initially in a af eigenstate, and the 
reservoir is in an arbitrary initial state. It turns out that 
the system’s behavior differs strongly between the resonance 

(|Ai?i—Ai? 2 | ^ I J|) and oMf-o/-reio«ance (I Ai?i—A£^ 2 1 3> 

|J|) regimes. To demonstrate this it is useful to transform 
to the rotating frame defined by \(j)) i—)■ exp[iA£^ 2 f(c’'f + 
erf)] \(j)). Both the leakage probability of the system qubit 
and the Heisenberg Hamiltonian are invariant under this uni¬ 
tary transformation. Thus, the new total Hamiltonian in 
the rotating frame can be obtained from iTtot(i) by replac¬ 
ing AEk AEk — AE 2 , k = 1,2 and /i 2 ,rest(i) '—>■ 
exp[iAi? 20 ’ 2 ]^ 2 ,rest(f) exp[—iAi? 20 'f]. Therefore, the sys¬ 
tem’s energy gap in this rotating frame is AEi — Ai? 2 - 

In the resonance regime leakage can occur in a time of 
0(1 J|“^), the fastest time allowed by the fundamental QSL 
bound Q. This happens, e.g., already in the case of single 
qubit reservoir, i.e., ft. 2 ,rest(i) = ^^rest(i) = 0, for which iTtot 
can easily be diagonalized. Under the resonance condition the 
states of the system and reservoir qubits are then swapped in 
a time of 0(| J|“^), so the fidelity with the initial state is lost. 
On the other hand, using our QSL bound (|^, we find that to 
have leakage with probability of 0(1) in the out-of-resonance 
regime, the minimum required time is lower bounded as 


^ I 71-1 n ~ ^ 

Tleak > c\j\ maxjl, ||, /,'|| } , 

maxt||/i2,i-est(u|| 


( 10 ) 


representing a potentially drastic increase in the time required 
for leakage relative to the minimum time c| J|“^ (where c is 
a constant of order one) obtained from more standard QSL 
bounds in the form of Eq. Q. 

This model has several interesting general implications: (i) 
Increasing the system gap can increase fidelity loss because 
the system may become resonant with reservoir DOFs; (ii) 
The relevant parameter which determines the speed of leak¬ 
age and fidelity loss is not the system gap but the energy 
mismatch between the system DOFs and the local reservoir 
DOFs, i.e., those that couple directly to the system. If they 
are in resonance with the system gap, then the reservoir can 
be insensitive to the gap, and leakage can happen in a time 
of order ~ ||Tfi|j“^, i.e., as fast as allowed by the fun¬ 
damental QSL bound (|^. On the other hand, if this en¬ 
ergy mismatch is large then relaxation is slow, even if the 
remote DOFs of the reservoir are in resonance with the sys¬ 
tem. (iii) Our QSL bounds are attainable in the regime where 


II [Tfi, iTR]|| ^ Ai?|jTfi||. (iv) Applying these bounds in dif¬ 
ferent rotating frames can lead to different independent con¬ 
straints. 

Beyond the AE > 2||iT/|| assumption. —^Finally, we dis¬ 
cuss how the large-gap assumption AE > 2||iTi||, used 
in deriving our previous bounds, can be relaxed. The key 
idea is to transform to a rotating frame in which AE be¬ 
comes larger. Let (Qq) be the projector onto the sub¬ 
space spanned by the eigenstates of Hs whose eigenvalues 
are greater (less) than those in I, and transform to the frame 
defined by |$) >->■ ) |<i)), where F S K. As we 

prove in the SM j^, the new gap between C and C-^ becomes 
AE + E. Consequently, bounds ([^-Q all remain valid for 
any E > 2||iTi|| — AE, after the substitutions 

AE^AE + F, H^it) ^ Hj,{t)-F{Q+ -Qc). (11) 

Moreover, using this generalization, we find that in the large 
F limit, bound Q implies that sin < 9||iTi|jf, which is 
the same as bound (|^, up to a constant. Thus, by varying F 
from 0 to 00 we can find a family of bounds which gradually 
changes from 0 to and find the strongest bound for fixed 
given values of the parameters. 

Conclusions .—In this work we introduced state- 
independent QSLs on leakage and fidelity loss in a 
Hamiltonian open system framework. The reservoir 
Hamiltonian iLR(f) only enters our bounds via iTR(f)]|j, 
implying that only local reservoir modes play a role in our 
QSLs. Another important conclusion concerns the common 
claim that increasing the system’s energy gap AE always 
results in better protection from coupling to the reservoir. 
The intuitive basis for this claim is the idea that a large 
gap suppresses thermal excitations by the Boltzmann factor 
g-AB/ZcT Under the BMA, the claim can be justified 
provided the spectral density of the reservoir is monotonically 
decreasing P^ . However, this condition is often violated, 
e.g., as in the case of an Ohmic bath. Our results, which 
are derived without approximations, demonstrate that this 
intuition is not always correct. Increasing AE can result in a 
resonance with the reservoir, causing the fidelity to drop on a 
timescale independent of AE, even if T = 0 and the reservoir 
is in a pure state. These results demonstrate the utility of 
state-independent QSL bounds for open system dynamics, 
and raise new questions about the efficacy of energy gaps in 
protecting quantum information. 
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Supplementary Material 


Appendix A: Closed system-Time dependent perturbations 


In the following theorem we hnd bounds on the effect of a time-dependent perturbation on the evolution of a closed system. 
This theorem is the main tool we use to prove bounds (0 and Q. The theorem has several other interesting applications which 
we shall explore in future work. The proof is given in two parts, in subsections |C2| and |C3| 

Theorem 1 Suppose the eigenvalues of the Hennitian operator Hq which are in the interval I C K are separated from the rest 
of the spectrum of Hq by at least AE. Suppose Pq the projector to the subspace spanned by the eigenstates of Hq whose 
corresponding eigenvalues are in T, and let Qq = I — Pq. Suppose V(t) is a time dependent perturbation that is turned on at 
t = 0 (i.e., V (t) = Ofor t < 0) and which for allt>0 is i) differentiable, ii) satisfies 2||l/(f)|j < AE. Let U{f) be the unitary 
evolution generated by the Hamiltonian H(t) = Hq + V(t), i.e., U(0) = I and dU(t)/dt = —iH{t)U{t). Then, 


\\QoUit)Po\\ < 2 


iino+)ii + iinoii 

AE 


/■'*' 2||l/(r)|| 

U AP-2|!V(r)|| ’ 


(Al) 


where || 1^ (O'*") || is the norm of perturbation at t = O'*". 

Furthermore, let 6E be the difference between the minimum and the maximum eigenvalues of Hq in X. Let Uoo (t) be the 
the unitary generated by the Hamiltonian H^[t) = PQH(t)PQ (or, equivalently, the unitary generated by H'^{t) = Hq + 
PQV{t)PQ), such that Uac,{0) = I, dUoc,{t)/dt = —iHao{t)Uao{t) (or dUaoif) / dt = —iH'^{f)Uao{t)). Then, the following 
bound holds 

||C/oo(t)Po - U{t)PQ\\ < 4 dr 

Note that the Hamiltonian H(^{t) = Hq + PqV (f)Po is the Hamiltonian obtained from first order perturbation theory in the 
limit AE —> oo. Bound ( |A2[ ) compares the dynamics generated by the actual Hamiltonian H{f) and the effective Hamiltonian 
dd'ooit)’ given that the initial state is in the support of Pq-, if the right-hand side is small then we know that for states which 
are initially in the support of Pq the final states under these two Hamiltonians are almost the same. The following Corollary 
highlights this point. 

Corollary 1 Let |'0) be an arbitrary state in the support of Pq. Let |'(/'oo(f)) = Uoo{t) |V') l'0(i')) = U{t) \'ip) be the states 

evolving under the evolutions generated by Hao(t) and H(t), respectively. Then the fidelity P(t) = \{'fit)\ipc>oit))\ between 
|'0oo(f)) <^nd \ fi(t)) satisfies the following bound: 


\\ViT)USE+\\ViT) 

AE 




AE-2\\V{t) 


(A2) 



^imo+)|| 

f\t\ 

/ rln- 

■||C(r)||(5P+||l/(r)||) ||l>(r)ll 

V2 

^ AE 

Jo+ 

AE AE-2\\V{t)\\ 


(A3) 


Proof. This Corollary follows directly from Theorem [T] together with the fact that 

< MW) - l^-W) II - III^W - 


< max I 


s/2 s/2 

\\[Uit)-U^{t)]\f)\\ _ \\[U{t)-U^it)]PQ\ 


|i#.)Gsupp(Po) ■ 


V2 


V2 


(A4) 


where in the first inequality we used -^/l — |((/)i|(/> 2 )| < (|| \4>i) — \(t> 2 ) ID/v/S, and in the last equality we used the definition of 
the operator norm. ■ 


Appendix B: From time-dependent closed systems to open systems (Proof of bounds (|^, (|^ and 

1. Proof idea 

To prove our bounds we consider the system and reservoir together as a closed system, and we transform to the rotating 
frame defined by the transformation i—> = ZR{t)\(j){t)), where Z^{t) denotes the unitary generated by the reservoir 
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Hamiltonian i7R(t), i.e. 


Z^{t) = iZ^{t)H^{t), Zr(0) = J. (B1) 

So, Z^it) = in the special case where H^(t) = H^, i.e., is time-independent. 

In this rotating frame the evolution of system and reservoir is governed by the Hamiltonian Hs -f where 

Hi{t) = Z^{t)HiZlit) . (B2) 

Note that the reduced state of the system does not change under this change of reference frame. In the rotating frame all the 
relevant information about the reservoir Hamiltonian and the interaction of the system with the reservoir is encoded in the 
term Thus, e.g., the difference between a reservoir with a “strong” Hamiltonian H^{t) and a reservoir with a “weak” 

Hamiltonian kH^{t) with A: ^ 1, is that in the case of the strong Hamiltonian Hi{t) oscillates much faster. Furthermore, any 
property of the reservoir Hamiltonian that does not show up in Hi{t) is irrelevant from the point of view of the system, i.e., the 
reduced state of the system remains independent of such a property. 

To prove bounds (|^ and Q we treat the interaction term Hi{t) as a time-dependent perturbation of the system Hamiltonian 
iFs- In particular, to prove bound Q we approximate the evolution generated by the Hamiltonian -f H\{t) with the evolution 
generated by the Hamiltonian obtained from first order perturbation theory, i.e., 

i7i(f) = Hs + (Pc 0 Ir) Hi{t) (Pc ® Ir) ■ (B3) 

Then we use Theorem to hnd an upper bound on the error we make because of this approximation. In particular, in the limit 
where the gap AE —> oo the approximation error vanishes. Indeed, with the exception of the term proportional to IS{PoHiPq), 
all the terms in bound Q arise from this approximation, and they appear because we use Hi (t) instead of the actual Hamiltonian 
Hs + Hi{t). The origin of the term involving IS(PoI^i^b) in bound Q is the second term in Hi, and it basically describes how 
the original evolution of the system, which is generated by Hs, is affected by the second term in Hi{t). In particular, if 
IS(PoI?i^o) = 0 then during the evolution generated by Hiit) the system is unaffected by the reservoir. 


2. Proof of bound 

Proof. First we move to the rotating frame dehned by the transformation 

o{t) d{t) = z^{t)o{t)zl^{t) 


(B4a) 

(B4b) 


and the hat denotes states and operators in the rotating frame. The system and reservoir together are a closed system and as a 
joint system they evolve unitarily. In the rotating frame |^(f)), the joint state of the system and reservoir, evolves according to 
the Schrodinger equation 


dt 


\i>{t))=-i{Hs + Hi{t))\^{i^) ■ 


(B5) 


The main idea is to treat Hi{t) as a time-dependent perturbation on Hs and to use the hrst part of Theorem[^ Let us write bound 
•lATJ as 


IIQo^W^oll < 2. 


ino+)ll + lint)ll 

AE 


2 ||l/(r)|| 


/o+ AE-2\\V{t)\\ 


(B6) 


where S{t) is the unitary generated by Hq + V{t), such that dS{t)/dt = —i{Ho + V{t))S{t) . 
To apply this bound we identify Hs ® /r as the initial Hamiltonian in Theorem[2 i.e.. 


Ho = Hs^ 


Pq = Pc Z) Ir 


and 


V{t) = Hi{t) . 


(B7) 


By dehnition S(t) is then the unitary generated by TJs ® -Ir + Hi (t), i-e., it is the unitary that describes the evolution of the system 
and reservoir in the rotating frame. This means that S{t) = Z^{t)U{t), where U{t) is the unitary describing the evolution of 
system in the lab frame, i.e., it is the solution of 17 = —i[Hs + H^{t) + Hi, U{t)], 17(0) = I. Then, using the fact that Hn{t) 
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commutes with Pc it follows that commutes with Pq, and so the left-hand side of bound (|B6[) is 


LHS = \\QoZ^{t)U{t)Po\\ = \\Z^{t)QoU{t)Po\\ = ||goC/(f)Po|| ■ 


(B8) 


On the other hand, using the facts that ||f?i(f)|| = ||iJi(f)|j and 
find that the right-hand side of bound (|B6[l is 


fAit) 


ZR{t)[Hj,{t),H,]zl{t) 


RHS 


4Hgi|| A , 2H[HR(f),gi]H 
AE AE-2\\Hi\\ 


||[i/R(f),i?i]||, we 


(B9) 


Therefore, bound (|B6|l implies 


\\QoU{t)Po\\ < 


Mm /■l‘l 2||[HR(f),gi] 

AE ^70+ AE-2\\H,\ 


(BIO) 


To prove bound (|^ we assume that the support of the initial reduced density operator of the system is restricted to C, i.e.. 


^0 0'sr(0)Po = {Pc Z> Ir) crsR(O) {Pc 0 Ir) = crsR(O) , 


(BID 


where (Tsr(O) is the initial joint system-reservoir state. Recall that the leakage probability is pieak(i) = Ti'[Qcf'(0]- B follows 
from bound ( |B10| ) that 


PieakD) = Tr [QoU{t)asR{Om{t)] = Tr [QoU{t)PoCTsR{0)PoUHt)Qo] 
< ||goD(DPol7lksR(0)||i 

^ f4\m , /■l*l , 2|i[iTR(D,i/i]||\' 

-^^ + D*AE-2||i/,|| ) ■ 


(B12a) 

(B12b) 

(B12c) 


where to obtain bound (|B12b|i we used the inequality |Tr(DB)| < HDjl ||B|| i, valid for any two operators A and B ll50l . ■ 


3. Proof of bound 0 

Proof. We shall use the second part of Theorem We consider the system and reservoir together as a closed system. Without 
loss of generality we assume that the state of the joint system is pure. If this is not the case then we can always purify the joint 
state by adding an auxiliary system which has a trivial Hamiltonian and is not interacting with the system and the reservoir. 

Let |$(0)) be the joint state of system and reservoir at f = 0, which need not be factorizable. Then, the joint state of system 
and reservoir at time t is given by |<h(f)) = U{t) |<1>(0)), where U{t) is the joint unitary evolution of the system and reservoir, 
and is generated by iJtot = Hs + E[R{t) + Hi, i.e. dU{t)/dt = —iHioi{t)U{t), and C/(0) = I. Thus |$(D) ^ purification of 

p{t). On the other hand, |4’(0)) is a purification of pii{t), i.e., 

TrR k(0)) ($(0)|e**^^) = pid(D . (B13) 


Next, consider the Hamiltonian 


Hi,,{t)=Hs + HR{t) + Pc®KR, (B14) 

where Kr G Herm('HR) is an arbitrary Hermitian operator acting on the reservoir Hilbert space, which will be determined later 
[below Eq. ( |B20[ )]. Note that this Hamiltonian can be obtained from the original Hamiltonian of the system and reservoir by 
replacing Hi with Pc Z) Kr. Also, note that from the point of view of states which are in the code subspace Pc, under evolution 
generated by this Hamiltonian the system and reservoir are decoupled from each other; Let C/dec(i) be the unitary evolution 
generated by Lfdec(i), such that dUdec{t)/dt = —iHdec{t)Udec{t), and C/dec(0) = I ■ Then, using the fact that Pc commutes with 
Hs, and H^{t) + Pc K^ acts trivially on |$(0)) we find 

TrR (c/dec(D |<i>(0)) {m\uL{t)) = TrR |<i>(0)) ($(0)|D‘"^) = p.d(t) (B15) 

In other words, Udec{t) |T’(0)) is another purification of pid{t). 
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Since U{t) |$(0)) and Udec{t) |‘&(0)) are, respectively, purifications of p{t) and pid{t), we can use Uhlmann’s theorem ll27l : 


F[pit),p,dit)] > ($( 0 )| Uljt)Uit) |$( 0 )) 


(B16) 


Define 


H^{t) = Hs + H^{t) + PoHiPo 

and let ?7oo(f) be the unitary generated by Hao{t), i.e. 

dUoo (f) 


Po = Pc ® Ir 


dt 


= -itHao{t)Uoo{t), t/oo(0) = I. 


Then, using the inequality < (II \4>i) ~ \4>2) ID/v^ and the triangle inequality we then find 


v/1 - F[p{t),p,d{t)] <^l- |($(0)| Ul^{t)U{t) |$(0)) 


< 


\\U{t) |$(0)) - Ud^t) |$(0)) 


^ \\U{t)Po-Ud,c{t)Po\\ ^ \\U{t)Po-U^{t) Po\\ 


V2 

-\\U^{t)Po-Ud,cit)Po\\ 


< 


V2 

\\U{t)Po - U^{t)Po\\ + ||C/oo(t) - C/dec(f)|| 

V2 


V2 


(B17) 

(B18) 

(B19a) 

(B19b) 

(B19c) 


where in inequality ( |B19b| l we used the definition of the operator norm (this is essentially identical to the proof of Corollary [TJ. 
Using Lemma [T] we have 


|Uoo(f) - Uda 


< 


(B20) 


f ds\\Hoo{s) - i7dec(s)|| = |f| \\PoHiPo -Pc® i^RlI • 

Jo+ 

Now we choose Kr to be the Hermitian operator for which \\PqHiPo — Pc ® -/TrU is minimized. Bound ( |B201 l thus implies 

WUooit) - C/dec(t)|| < |f| mmK,^n.Mn.)\\PoHiPo -Pc® Kr\\ < \t\ lS{PoHiPo) , (B21) 

where we used definition •H). As we mentioned in the main text, note that IS{PoHiPq) can be nonzero only when C has 
dimension larger than one. To see this note that if dimC = 1 then Pc = IV'XV'I for a normalized state \ip)-, then, choosing 
Kr = {tj;\Hi\tj;) we find (|X’)(V^| ® Ir)F[i{\iP){ ip\ ® Ir) = |■*/')(V'I ® (X’l'ffilX’)^ so IS(Pof7i^o) = 0 in this case. 


Using bound ( |B21| ) in bound ( |B19| l, we find 


Vl-F[p{t),p,dit)] < ^oo(f)Po|| ^ ^isiPoHiPo) . 


We show below, using bound ( |A2| ) of Theorem[T] that 

4||f7i|| , 


\\U{t)Po-U^it)Po\\ < 


AE 




/■l*l ||[gi,iJR(r)]|| 
U A£;-2||iT,|| 


(B22) 


(B23) 


Combining this with bound ( |B22| l proves bound Q. It thus remains to prove bound ( |B23| l, which we do next. 

The argument is analogous to the argument we used to prove bound (|^ using bound (|A1|). As in Sec. B2 


we start by 


transforming to the rotating frame defined by |'I'(f)) i—)• |4'(f)) = ^R(f) |^'(f)), where ZR(t) is defined in Eq. (BTJi. In this 
rotating frame |'P(f)) (the joint state of the system and reservoir) evolves according to the Schrodinger equation l'I'(f)) = 

—i{Hs + I'f (f)), where Hi{t) = ZR{t)H\Zldt). Again, the idea is to treat H\{t) as a time-dependent perturbation on iJs 

and to use bound ( |A2| i of Theorem[2to estimate the effect of this perturbation. Let us rewrite bound ( |A2| i in the form 


\\S^{t)Pd-S{t)Pd\\ < 


4||V(0+)|| 

AE 




f.|t| n\t\ 

dT\\V{T)\\{6E+\\V{T)\\)+A dr 


l|f>(T)|| 


b+ 


0+ AL;-2||U(r)r 


(B24) 


where S{t) and Sao{t) are, respectively, the unitaries generated by the Hamiltonians Hq + V (f) and Hq + PqV{ t)Po. 
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To apply bound ( |A2[ ) we identify Hs 0 Ir as the initial Hamiltonian Hq in Theorem[2 i.e., 

Ho = Hs®Ir, Po = ^C0/r and V{t)=Hi{t). (B25) 

First, consider the left-hand side of Eq. ( |B24| i. Since S{t) is the unitary generated by -f Hi{t), then S{t) = ZR{t)U{t). 
Similarly, Soo {t) is the unitary generated by 

Hs + (Pc 0 lR)Hi{t){Pc 0 Jr) =Hs + ZR{t) (PoHiPo) 4(f) . (B26) 


Going from the rotating frame to the lab frame, one can easily see that this Hamiltonian transforms to Hs + Hj^ + PqHiPq 
which, by definition [Eq. ( |B17| l] is Hao- Since Uoo{t) is the unitary generated by Hao, it follows that Soo{t) = ZR(t)Uoo{t)- 
Therefore, the left-hand side of Eq. ( |B24| l is 

LHS = \\S^it)Po - ^(t)Poll = II^R(f) [U{t)Po - U^{t)Po] II = \\Uit)Po - U^it)Po\\ . (Bll) 

On the other hand, using Eq. (|B25[l we can easily see that the right-hand side of Eq. (|B24|l is 


RHS 


4||gi|| 

AE 


4 

AS 



dT\\Hi\\{SE + \\Hi\\) + 4 



\\[H,,Hr{t)]\ 

AE-2\\Hj\\ 


4\\H,\\ , ^^JH,\\{SE+\\H,\\) , ^ ||[S,,iTR(r)]|| 

-AE-+ ‘‘S ‘*Ae-2||H,|| 


(B28) 

(B29) 


where again we used 114(^)11 = ||-ffi||, and \\-^Hi{t)\\ = ||[Sj, iTR(f)]||. 
Substituting Eqs. ( |B27| i and ( |B29[ ) into bound ( |B24| i we find 


||C/(t)Po-C/oo(f)-Po|| < 


4||gi|| , ^,^, ||giH(^S+||gi||) 
AS ' ' AS 



||[Si,SR(r)]|| 

AS- 2 ||Si|| 


which proves Eq. (|B23|l, and so completes the proof Eq. Q. ■ 


4. Proof of bound 


Proof. We present the proof for the special case where HR{t) = Hr is time-independent. The time-dependent case follows 
exactly in the same fashion. 

Let |$(0)) be any arbitrary joint pure state of the system and reservoir at f = 0. Then, the joint state of the system and 
reservoir at time t is given by |$(f)) = |$(0)), where Stot = Hs + Hr -f Hi. So, by definition, |'i>(f)) is a purification 

of p{t), the reduced state of system at time t, i.e.. 


p{t) = TrR{mt)) mm- 


On the other hand, one can easily see that |4>id(f)) = e |(j)(o)) is a purification of pid{t), i.e., 

TrR (|<l>.d(f)) ($.d(f)|) = e-*‘^^TrR (|$(0)) ($(0)|) = p.^t) . 

Therefore, using Uhlmann’s theorem we find 


F[p(f),Ad(f)]>|(‘i>(t)|$M(f))|. 


Next, using the fact that \J\ — |((/)i \(f)2) \ < (|| |^i) — e*® 
E[p{t),pid{t)\ < 4l - \{4>{t) |$id(f))| < ■ 

_ g-2Ctg 

“■ TT 


im \\)/V 2 , where 6 is an arbitrary phase, we find 

|||$(f))-e—*| 4 >.d( 0 )|| 

72 


-it(Hs+H^) |$(o)) 11 


||g-itrr,o, _ g-ictg-it(rrs-i-ffR) |j 

' 71 


(B30) 

(B31) 

(B32) 

(B33a) 

(B33b) 



































11 


where c is an arbitrary real constant. Finally, using Lemma[T]we have 

_g-ictg-zt(rrs+i/R)|| < |i| \\H^^_^Hs+HT, + cI)\\ = \t\ \\Hi-cI\\ < |t|||iF,|| , (B34) 

Thus, \/l — F[p{t), pid{t)] < 1^1 \\Hi — cI\\/\/2 for any c S K. Let Amax and Amin be the maximum and minimum eigenvalues 
of Hi respectively. Thus the maximum and minimum eigenvalues of Hi — cl are Amax — c and Amin — c, respectively. It is easy 
to see that to minimize ||7Ji — cl\\ we should choose c = (Amax + Amin)/2, whence ||iFi — cl\\ = (Amax — Amin)/2. Therefore 
y/l - F[p{t),pid{t)]/y/2 < |f|(Amax “ Amin)/4 < |f11|iFi|1/2, as claimed. ■ 


5. Proof of transformation i fTT) 

Recall that I is an interval of K which includes at least one eigenvalue of Hs, C is the subspace spanned by the eigenvectors 
of Hs whose eigenvalues are in X, and Pc is the projector onto the subspace C. Recall that AE denotes the gap between the 
energy levels of Hs inside and outside C (i.e., if Ai and A 2 are two distinct eigenvalues of Hs such that Ai G I but A 2 ^ X, then 
|Ai - A 2 I > AE). 

Proof. Let Qc (Qc ) be the projector onto the subspace spanned by the eigenstates of Hs whose eigenvalues are greater (less) 
than those in X. This definition implies that 


Q^ + Qc+Pc=h- (B35) 

To prove bounds (|^-Q, we moved to the rotating frame described by the transformation 

|T'(f)) ^ |^(f)) = Zr(<)|^(<)), (B36) 

where Zii{t) is defined in Eq. ( |B1| |. To prove transformation ( [TT] i we move to a new rotating frame defined by 

|T'(f)) ^ |T'(f)) = Wpit) |T-(f)), (B37a) 

WF{t) = exp{-itE[Q+ - Qc ]) ® (B37b) 

(i.e.. Writ) is generated by Hw{t) = H^{t) — E[Q'^ — ]), where F is an arbitrary real number satisfying 

L’> 2||iTi|| - at; . (B38) 


We will shortly present the motivation for this condition (more generally, F can be chosen to be time-dependent, but we shall not 
consider this here). In the rotating frame |'L(t)) (the joint state of the system and reservoir) evolves according to the Schrodinger 
equation 


^^\^,(t))=-^{Hll + Hl{t))\^{t)) , 

where 


(B39) 


Hii = Hs+ F[Q+ - Qc] and Hi{t) = WF(t)HiWl{t) . 


(B40) 


Thus |^(f)) = S{t) l'I'(O)), where S{t) is the unitary generated by Hq + Hi{t). Note that S{t) can be written as S{t) = 
WF{t)U{t), where U{t) is the unitary describing the joint evolution of the system and the reservoir in the lab frame, i.e., the 
solution of dU{t)/dt = —i[Hs + Hfi{t) + Hi]U{t) with U{0) = L 
Again, we use the first part of Theorem [T] Let us rewrite bound ( |A1| | as 


\\QoS{t)Po\\ < 2 


Hgi(0+)|| + ||gi(f)|| 

AD 


+ 


Jo+ AD - 2\\HiiT)\\ 


(B41) 


where Pq = Pq ® /r. Let AD denote the gap of Hamiltonian Hq. The conditions of Theoremj^require that AD > 2||-ffi(r)|| 
for all 0 < T < t, and we next show that this can be satisfied provided we choose F as in condition ( |B381 l. 

We claim that the gap AD of Hq is AD = AE + F. To see this, first note that by definition and Q'^ are diagonal 
in the eigenbasis of Hs, so that Hq is diagonal in the same basis. Adding FlQ'^ — (5(7] to Hs has the effect of adding F to 
all the eigenvalues of Hs greater than those in X, and subtracting F from all the eigenvalues of Hs less than those in X, while 
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leaving the eigenvalues in I alone. Since the gap between the eigenvalues in X and those not in I was AE' before the addition of 
F[Q'^ — Qq], it becomes AE + E after this addition, which is then the new gap between C and the orthogonal subspace. Then, 
since ||E^[(<)|| = ||Ei|j, condition ( |B38| l guarantees that 

AE = E + AE>2||Ei|| , (B42) 


and so we can apply Theorem [T] and bound ( |B41| i. 


Using the fact that Hw{t) = H^{t) — F[Q'^ — ] commutes with Pc it follows that Wpit) commutes with Pq, and so the 

left-hand side of Eq. ( |B41| l is 


LHS = IIQo^W^^oll = \\QoWF{t)U{t)P4 = ||QoC/(i)-Po|| 

On the other hand, using the facts that ||Ei(f)|| = ||Ei|| and 

= jwF{t)[Hi,Hwit)]wl{t)j = ||[Ei,ER(f)-E(Q+-Qc)]| 
we find that the right hand side of bound ( |B41| l is 

4||E,|| , fW ^^2\\[H^,H^it)-FiQ+-Qc)]\\ 


RHS = 


AE-fE 


/o+ 


AE-f E- 2||Ei| 


Therefore, bound (|B41|i implies 


\QoU{t)Po\\ < 


Mm\ , /■l‘l^^2||[Ei,ER(f)-E(Q+-Qc)]| 


AE-f E 


/o+ 


AE + F-2\\Hi\' 


(B43) 


(B44) 


(B45) 


(B46) 


Suppose (Tsr(O) is the initial joint state of the system and reservoir with the property that the reduced state of the system has 
support only in C, i.e.. 


Pq o'sr(0)Po — {Pc ® -^r) ctsr(O) {Pc ® Ir) — o’sr(O) ■ 

Then, bound ( |B46| ) implies 

Pieak(f) = Tr [C/(f)<TsR(0)E'l'(f)Qo] 

= Tr [QoU{t)PoasR{0)PoU\t)Qo] 

< IIQoC/(i)i"of lksR(0)||i 

^ / 4||Ei|| /-'‘I ^ j\\[Hi,H^{t)-F{Q+-Qc)] 

-yAE + F^Jo^ ^ AE + E-2||Ei|| 

which is bound (|^ after the replacements 

H^{t) ^ H^it) - F{Q+- Qc) 

AE 1-^ AE -f E . 


(B47) 


(B48a) 

(B48b) 

(B48c) 

(B48d) 


The argument which leads to bound Q can be repeated in a similar fashion. ■ 


(B49a) 

(B49b) 
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Appendix C: Proof of theoremj^ 
1. Preliminaries 


a. Three useful Lemmas 


Lemma 1 Suppose the unitaries Ui{t) and U 2 {t) are generated by the Hamiltonians Hi{t) and H 2 {t), such that 17i,2(0) = I 
and dU 1 ^ 2 {t)/dr = — Then, 


f\t\ 

m{t) - U 2 m < dr - H 2 iT)\\ . (Cl) 

Jo 

See Sec. |D l| for the proof. 

Lemma 2 Let P{t) be a time-dependent projector satisfying [P{t), H(t)] = 0 for all times t, where H{f) is the Hamiltonian 
generating the unitary U(t), i.e., U{t) = —iH{t)U(t) with (7(0) = I. Let H,runc{t) = P{t)H(t)P{t) be the generator of 
Ufruncit), be., dU truncif) / dt = —iHtnmc{t)Utrunc{t) with UtrunciO) = I ■ Assume the time-derivative of P{t) exists for ( > 0 . 
Then 


\\U{t)P{Q+)-P[t)U{t)\\ < / ||P(r)||dr 


/o+ 

\\U{t)P{0+) - Urrur,c{t)P{0+)\\ < 2 f ||P(r)||dr 

7o+ 


See Sec. D 2 for the proof. 

Lemma 3 Let P and P be the projectors with the same rank. Then 

\\P - P\\ = \\PQ\\ = \\PQ\\ 


(C2a) 

(C2b) 


(C3) 


where Q = I — P and Q = I — P. 


See Sec. D 3 for the proof 


b. Bounds on the effect of perturbations 


Let Iq ^ be interval containing one or more eigenvalues of a Hermitian operator Hq. Suppose the eigenvalues of Hq 
in the interval Iq are separated from the rest of the eigenvalues of Hq by at least AE. I.e., for any pair of eigenvalues Ai and A 2 
of Hq, if Ai G Iq and A 2 ^ Iq, then |Ai — A 2 I > AE. 

Lemma 4 Let V be a Hermitian operator satisfying 0 < ||17|| < AE 12. Let I be the interval obtained from Iq by adding 
a margin of ||lf || on the left and the right of Iq. Let Pq and P be, respectively, the projectors onto the subspaces spanned by 
eigenstates of Hq in the interval Iq and eigenstates of Hq + V in the interval I. Then, Pq and P have the same rank and 


\\P 


^oll < 


2M 

AE 


(C4) 


This is Lemma 3.1 of Ref. ifSTl . (See also Theorem VII.3.2 of Ref. ifSOl ). 

Lemma 5 LetV{f) be an arbitrary differentiable Hermitian operator satisfying Q < ||C(t)|| < AE 12. Letl{t) be the interval 
obtained from Iq by adding a margin 0 / ||C(()|| on the left and the right of Iq. Let P{t) be the projector onto the subspace 
spanned by the eigenvectors of Hq + V (t) whose eigenvalues belong to I{t). Then, 


^ 2\\P{t)V{t)Q{t)\\ 

A£;-2||V(()|| 

(C5a) 

/ 2|iV(()|| 

- AE-2\\V{t)\\ ■ 

(C5b) 


See Sec. D 4 for the proof 
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2. Proof of the first part of Theoremj^ [bound jAlU 


First, using the triangle inequality, we note that for an arbitrary projector A it holds that 

||(/ - A)17(t)A|| < ||(/ - A)t/(t)P(0+)|| + 11(7 - A)U{t) [A - P(0+)] || (C6a) 

< 11(7 - A)U{t)P{0+)\\ + ||A - P(0+)|| , (C6b) 

where to get the last inequality we used the fact that 7 — A is also a projector and so ||7 — A|| = 1, together with the unitary 
invariance of the operator norm. Similarly, we find 

||(7-A)P(f)P(0+)|| < ||[7-P(f)]P(f)P(0+)|i + ||[P(f)-A]P(f)P(0+)|! (C7a) 

< ||[7-P(f)]P(f)P(0+)|! + |lP(f)-A|i. (C7b) 

Combining these two bounds we find 

11(7 - A)P(t)A|| < ||[7 - P(f)]P(f)P(0+)|| + ||P(f) - All + ||A - P(0+)|| . (C8a) 

Then, we observe that 

||[7-P(f)]P(f)P(0+)|| = ||P(f)P(0+)-P(t)P(t)P(0+)|| = 11 [P(f)P(0+)-P(f)P(t)]P(0+)|| (C9a) 

<\\U{t)P{0+)-P{t)U{t)\\. (C9b) 

Substituting this into bound ( |C8a| ) we find 

11(7 - A)P(f)A|| < ||P(f)P(0+) - pmm + \\P{t) - All + ||A - P(0+)|1 . (ClOa) 


The above bound holds for any projector A. Choosing A = Pq, we can easily see that bound ( |A1| | follows immediately using 
this bound together with lemmas and 


3. Proof of the second part of Theoremj^ [bound i |A2[ (] 


Proof. First, without loss of generality we assume the minimum energy in 2 is zero, and so the maximum eigenvalue in 2 is 
(5P. ll52l Let 2(t) C K be the interval obtained from 2 by adding the margin ||L(f)|| on the left and the right of 2. Let P(t) be 
the projector onto the subspace spanned by the eigenstates of H{t) whose eigenvalues are in 2{t). Then, since by assumption 
2|| L(f) II < AE, it follows from Lemma|^that the rank of P{t) is the same as the rank of Pq, the projector onto the eigenstates 
of 77o with eigenvalues inP. Furthermore, Lemmas and imply 

||Q(f)Po|| = ||P(f)Qol! = l|i^(t) - J^oll < , (Cll) 

where Q{t) = I — P{t) and Qq = I — Pq . 

Define the Hamiltonian 77tnjnc(f) to be the truncated version of H{t) in which we have removed all the energies of 77(f) 
outside 2{t), i.e.. 


77trunc(i) = P{t)H{t)P{t) . (C12) 

Let Ptrunc(f) be the unitary generated by 77t,unc(f), i-e., dUtmadt)/dt = -*77trunc(i)C7trunc(i), and Ptrunc(O) = 7. Let P(0+) be 
the projector onto the subspace spanned by the eigenstates of 77(0’*') = Hq + 1^(0+) whose eigenvalues are in P(0+), where 
Hq + V (0'*') is the Hamiltonian of the system immediately after the perturbation is turned on. Then, we have 

||Loo(f)Po - L(f)Po|| < ||Poo(f)Po - Utmndt)Po\\ + ||P,™nc(i)Po ” C/(t)Po|l 

< ||C/oo(f) - t/,^„c(t)|| + II [U{t) - U^,dt)] [^(0+) + {Po - P(0+))] Poll 

< ||C/oo(f) - 7/^nc(t)|| + II [U{t) - Ptrunc(t)] P(0+)|| + 2||Po - P(0+)|| . 


(Cl 3a) 
(C13b) 
(Cl 3c) 
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Using Eq. ( |C1 l| l we find 

\\U^{t)Po - U{t)Po\\ < ||t/oo(f) - U,^„c(f)|| + \\U{t)P{0+) - C/t„„c(t)P(0+)|| + . (C14) 

Lemma [^already provides us with an upper bound on ||f/(f)P(0“*') — C/trunc(f)^’(0^)||. Thus, to prove bound ( |A2| l it remains to 
hnd an upper bound on ||?7trunc(f) — C4o(f)||- Indeed, we prove below that 


||I7tru„c(f) - C/oo(f)|| 



||U(r)|| {6E+\\V{r)\\) 
AE 


Substituting this bound and bound (|C2b|) into inequality (|C14|l we hnd 


(C15) 


\\U^{t)Po-U{t)Po\\ 


- 1 /'"ir Hn^)ll(^-5;+||U(r)H) 
■ 7o+ al; 



dT\\P{T)\\ + 


4||l"(0+)li 

AE 


(C16) 


Combining this with Lemma]^ which puts abound on the norm of ||P(T)|j, proves the second part of Theorem[l] It thus remains 
to prove bound ( |C15| l, which we do next. 

To prove bound ( |C15| l, we recall that Uoo{t) and C/trunc(f) are the unitaries generated by PQE[{t)Po and P{t)H{t)P{t) respec¬ 
tively. So, we hrst hnd a bound on the difference of these two Hamiltonians. To do this, we note that 


P{t)H{t)P{t) = P(t)H{t) = P{t)H{t)Po + P{t)H{t)Qo = P{t)H{t)Po + P{t)H{t)P{t)Qo . (C17) 

By dehnition P{t) is the projector onto the subspace spanned by the eigenstates of H{t) whose eigenvalues are in I{t). Recall 
that 2{t) is the interval obtained by adding the margin of || U(f) || to the interval I. So, it follows that all the eigenvalues of El (f) 
in 1(f) are between — ||U(f)|| and SE + ||U(f)||. Thus all the eigenvalues of P{t)H{t)P{t) are likewise between — ||U(f)|| and 
6E + |jU(t)||, and so \\P{t)E[{t)P{t)\\ < 6E + ||U(f)||. Therefore, using Eq. ( |C17| i we hnd 



\\P{t)H{t)P{t)-P{t)H{t)Po\\ = \\Pm{t)P{t)Qo\\ < \\Pm{t)P{t)\\\\Pit)Qo\\ 

(Cl 8a) 


<{SE+\\Vm)\\Pit)Qo\\- 

(Cl 8b) 

Next, we observe that 



PoHit)Po = Pit)H{t)Po + [Po - Pit)]H{t)Po . 

(C19) 

Therefore 

\\PoH{t)Po - P{t)H{t)Po\\ = II [Po - Pit)] [U(f)Po + HoPo] || 

(C20a) 


< ||Po-P(f)||||U(f)Po + PoPo|| 

(C20b) 


< ||Po-P(t)||(||U(f)|| + ||PoPoll) 

(C20c) 


< ||Po-P(t)||(||U(f)||+^P) 

(C20d) 


= \\Pit)Qo\\i\\V{t)\\+6E), 

(C20e) 


where to get inequality ( |C20d| l we used the fact that Pq is the projector onto the eigenstates of Hq with energy in Xq, which is 
between 0 and SE, and to get inequality ( |C20e| l we used Lemma 
Combining bounds ( |C18| ) and ( |C20| l we hnd 

\\PoH{t)Po - Pit)Hit)P{t)\\ < \\PoH{t)Po - Pit)H{t)Po\\ + \\P{t)Hit)P{t) - P(f)i7(f)Po|| (C21a) 

<2{\\V{t)\\+6E)\\P{t)Qo\\ . (C21b) 


Using Eq. ( |C1 l| l we have ||P(f)Qoll < therefore 


\\PoH{t)Po-P{t)H{t)P{t)\\ 


^ ||u(f)||(JL; + ||u(f)||) 
AE 


(C22) 


Einally, since P{t)E[{t)P{t) generates the unitary Utrunc(f) and PoH{t)Po generates the unitary Uoo{t) with the initial condi¬ 
tions Uoo(O) = Utrunc(O) = I, using Lemmaj^we have \\U^{t)Po - U(f)Poll < /d*' ll^o^^('r)Tb - PiT)H{T)P{T)\\, and 
bound (|C15|) follows, as claimed. ■ 
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Appendix D: Proof of the preliminary lemmas 
1. Proof of Lemmajl] 


The proof of this Lemma is given, e.g., in Ref ll5^ and is reproduced here for completeness. 

Proof. The proof is a straightforward application of the unitary invariance of the operator norm together with the triangle 
inequality: 


\mt)-u2m = \\ukt)ui{t)-i\\ 

< ||^‘rfs^(C/|(s)C/i(s))|| 


< 



Ul{s)H2{s)Ui{s) - ul{s)H,{s)Ui{s) || 


< 


< 



uUs)H2{s)Ui{s) 



ds\\H2{s) 




UUs)H,is)U^{s) II 


(Dla) 

(Dlb) 

(Die) 

(Did) 

(Die) 


2. Proof of LemmaH] 


a. Proof of bound \C2&) 


Proof. Let S{t) = W{t)P{t)U{t). Then 

\\U{t)P{0+) - pmm = \\UHt)Pm{t) - P(0+)|| = ll^(t) - 5(0+)|| = f S{t) dr 

J 0+ 

(iU^^( t)F{t)P{t)U{t) (t)PU{t) - iU\T)P{T)H{T)U{T)^ dr 
U\t)P{t)U{t) dr 


10+ 

pt 


10+ 

.|t| 


< 


/o+ 


\\UHt)P{t)U{t)\\ dT= f ||P(t)|| dr 
Jo+ 


where to get Eq. (|D2c|l we used the premise that [P{t),H{t)] = 0 Vt. 


(D2a) 

(D2b) 

(D2c) 

(D2d) 


b. Proof of bound l |C^ 


Proof. First, using the triangle inequality we find 


||C/(f)P(0+) - C/t.-unc(t)P(0+)|| < ||17(t)P(0+) - P(f)t/,runc(t)|| + || P(t) C/trunc (t) " (f)P(0+) || 

= \\U\t)P{t)U,mnc{t)- P{0+)\\+\\ 

^trunc (t)P(t)C/,„„c(f)-P(0+)|| . 


By bound (|C2a|l we already know that the second term satisfies 


|t| 


||P(r)||dr . 


(D3a) 

(D3b) 


11 ^trunc (f)P(f)17„unc(i)-P(0+)|| < 


0 + 


(D4) 
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Now let S{t) = {t)P{t)Uauncit)- Then 


\\uHt)p{tn,,,,{t) - P(0+)|| = \\s{t) - 5(0+)|| = 

Ht| 


n|t| 


10 + 


S{t) dr 


(D5a) 


(^C/t(r)i^(r)P(T)C7.„„e(T) + UHt)P{t)U^,,{t) - ^C/t(r)P(T)i^.„„e(T)C/tm„c(r)) 

(D5b) 

^zC/'l'(r)iJtrunc('r)17ti-unc(T) + 17''’(T)P(r) Ptrunc (t) - zC/'''(T)Ptrunc (l)17ti-unc (t)) dr 

rm . f\t\ 

< \\U^T)P{T)Utrunc{T)\\ dr = \\P{T)\\dT, 

Jo+ Jo+ 


' 0 + 

r\t\ 

lo+ 


(D5c) 

(D5d) 


where to get bound ( |D5c| i we used the definition Ptrunc('r) = P{t)H{t)P{t) = P{t)H(t). Combining these bounds we have the 
claimed result: 


||P(i)P(0+)-Ptrunc(l)P(0+)| 


r\t\ 

<2 \\P{r)\\dr. 

Jo+ 


(D6) 


3. Proof of Lemma|3] 


Proof. First note that since the supports of P and P have the same dimension, then there exists a unitary U such that UPU^ = P, 
and therefore UQW — Q. This implies 

||QP|| = ||QPPpt|| = ||Q[/p|j. (D7) 

Using the fact that P = PW (P + Q)UP we find that 

IIP - PU^PUPW = \\PU^QUP\\ = IIQUPf = IIQPPP^f = IIQPf , (D8) 

where we have used the fact that ||AAl|| = || A|p for any operator A. Next, we use the fact that for any operator A, the operators 
AA'' and A'''A have the same eigenvalues. This implies that PU^PUP and PUPWP have the same eigenvalues. It follows 
that 


IIP - PPPP^PII = IIP - PP^PPPII = IIQPf. (D9) 

Then, using the fact that P = PU{P + Q)WP we find 

||ppgp''p|| = IIP - PPPP^PII = IIQPf, (DIO) 

which implies 

\\PUQ\\ = ^fipUQU^\ = \\QP\\. (Dll) 

The left-hand side is equal to ||PPQ|| = ||PPQP'1|| = ||PQ||. Therefore, we find ||PQ|| = ||QF’||. 

To prove that IIP — P|| = ||PQ|| = ||QP||, note that 

IIP - P|| = IIP _ (P + Q)P|| = IIP _ PP _ QP|| = ||PQ _ gp|| = ^\\QPQ + PQP\\ , (D12a) 

where to get the last equality we have used the fact that for any operator A, \\A\\ = Since the supports of PQP and 


dr 
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QPQ are orthogonal, it follows that 

yJ\\QPQ + PQP\\ = ^Jm?ix{\\QPQl\\PQP\\} = m?ix{^J\\QPQl \J\\PQP\\} = max{|lPQ||, ||PQ|1} . (D13) 

Since ||P(5|| = H^’QH we find that 

IIP - P|| = max{||PQ||, IIPQII} = ||PQ|| = ||PQ|1 . (D14) 


4. Proof of LemmaH] 


Proof. Recall that for any pair of eigenvalues Ai and A 2 of Hq, if Ai G Pq ^2 ^ Po> then |Ai — A 2 I > AP. Therefore, since 
||lA(t)|| < AP/2, the number of orthonormal eigenstates of Hq + V (f) in I is equal to the number of orthonormal eigenstates 
of Hq in Iq. Thus the rank of P(t) is time-independent. Furthermore, the eigenvalues of Hq + V (t) in I are separated from the 
rest of the spectrum of Hq + V (t) by at least AE — 2||lA(f)||. Let Ptot(f) = Hq + V{t), APtot(f) = Hiot{t + At) — Ptot(f) and 


Thus 


= P(f)APtot(f)P(f) + Q{t)AH,,,{t)Q{t) (D15) 

AH^°^\t) = P(f)AP,ot(i)Q(f) + g(f)AP,ot(f)-P(i) (D16) 


+ At) = 




it). 


(D17) 


The terms inside the bracket are block-diagonal with respect to P{t) and Q{t). So, in the limit where At is sufficiently small 
such that II AP(o°‘'‘®^(f)|| < AE — 2||14(f)||, we find that P(t) is also the projector onto the subspace spanned by the eigenstates 
of Ptot(i) + (f) . Note that the eigenvalues of Ptot(i) + ^Hi^'’^^\t) whose eigenvectors are inside the support of 

P{t) and the eigenvalues whose eigenvectors are outside the support of P{t) are separated from each other by at least AE — 
2||y(f)||-||Ap(°“®)(f^|. 

Next, we use Lemma|^to find the effect of adding to H^otit) + According to Lemma|^ 

||P(i + A,)-PW||< _2||AHrWII , 2||PmAH,.,WQW|| 

AB - 2||r(()|| - II Aff,5™>(()ll AE - 2||r(()|| - || AB<“''(t)|| 

In the limit where A goes to zero, this implies 

„ ^ 2||P(f)P.o,(t)Q(f)|| _ 2||P(f)l/(f)Q(f)|l 

" ^ - AP-2||f^(f)|| AE-2\\Vit)\\ ■ ^ ’ 

This proves bound ( |C5a| i. Bound ( |C5b| i follows since the largest eigenvalue of the off-diagonal part of ||1A|| cannot be larger than 
the largest eigenvalue of |jl/||. 

We note that, alternatively, one can prove Lemma|^using standard resolvent techniques, i.e., using P{t) = - 2 ^ /r P(z, t)dz, 
and P{t) = — 5 ^ fp H(z, t)V(t)R(z, t)dz, where Riz, t) = (Hq + V{t) — zl)~^ is the resolvent, and T is a properly chosen 
contour in the complex plane around X{t). See, e.g., Il54l . ■ 


Appendix E: Proof of remaining statements in the main text 
1. The rate of energy exchange between the system and the reservoir and the interpretation Q,{t) 

In the main text we claimed that in the special case where 

1. H^(t) = Pr and hence also f2(f) = 12 (time-independent), and 

2. C is the bottom (top) energy sector, i.e., all energy levels in C-^ have energy less (larger) than energy levels in C, 
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we can interpret 17“^ as the minimum time the reservoir needs to transfer the required energy to move the system from C to C^. 
Here we present the argument in more details. 

The key point in the argument is that in the time-independent case the total energy of system and reservoir together is a 
conserved quantity. So, if the energy of the system changes considerably, the energy of reservoir {Eb) = Tr[/5sR(f)^^R] should 
also change, where psR(f) is the joint state. But the maximum rate of change of energy of reservoir is equal |^(f5B)| < 
II [iJi, Hr\\\, with equality for some psR{t). (This follows easily from multiplying psR = — p] by taking the trace and 
using its cyclic property to cancel the terms not involving Hi, and using the inequality |Tr(yli3)| < ||A|| ||i3||i ISOl .) 

Next we notice that the second assumption implies that to move the system state from C to C^, the reservoir energy should 
change at least by IS.E — 2||i?i|| (note that adding Hi to Hs can shift each eigenvalue of iJs by at most ||Tfi||, and so the gap can 
shrink to /S.E — 2||iJi|j). Since ||[iTi, Hn]\\ is the maximum rate of change of {Eb), exchanging this amount of energy takes a 
time of at least (AE’ — 2||iJi||)/|| \Hi, iTR]||. So, under these assumptions 17“^ (up to a factor of 1/2) can be interpreted as the 
minimum time it takes the reservoir to transfer (absorb) the required energy to go from C to C^. This explains the factor fH in 
bound (|^ [and bound Q] in the time-independent case. 

However, this simple argument for the minimum energy exchange time cannot explain some important aspects of bounds (|^ 
and Q. In particular, it does not apply when C is neither the bottom nor the top energy sector. In this case the leakage can 
happen without any change in the average energy of the system and so one might expect that leakage can happen in a much 
shorter time. However, bound (|^ shows that this is not the case. Nor does the simple argument explain the fact that the bound 
on the probability is quadratic in ||Tfi||, or the fact that the bound holds even for time-dependent Hamiltonians, where the energy 
is not generally conserved. 


2. Modification of by the reservoir as a beneficial effect 


In the main text we commented that in some cases the modification of H^ by the reservoir can be a beneficial effect. This 
happens when PqHiPq = PcK^Pc ® Ir for some Hermitian system operator Ks. In this case, in the AE oo limit the 
evolution of states inside C is described by the Hamiltonian Hs -f Ks. This modification of the system Hamiltonian can be 
useful, e.g., in the context of quantum computation driven by dissipation ll43l . We can use the argument leading to the bound 
to find how close the actual evolution is to the ideal evolution generated by Hs -I- Ks- Let pid{t) = 

Then, it turns out that for any state p(0) G C: 


a/1 - E{pid{t),p{t)) < 


A M 

A [ AE 


H[gl,gR]H 

AE-2||Ei|| 


\Hi\\{6E+\\Hi\ 

AE 


(El) 


where we considered the case of time-independent Er for simplicity. We note that the term appear in the 

Markovian analysis of Ref. Il43l . 


3 . Proof of the inequality E[p(i),pid(f)] < ^/K^p]edk(t) 

Recall that we defined the leakage probability as pieak(f) = T'^[p{t)Qc]- Consider two states cr [i.e., pid{t)] and t [i.e., p{t)] 
where PccrPc = cr for some projector Pq {Pc = Qc = I ~ Pc) with support C. Considering the Taylor expansion of ^/(T, 
we can easily see that 

= ^/^Pc = Pcy/a ■ (E2) 

Define r' = PctPc/Tt{Pct). Then, 

E(r, cr) = Trj^ = Trj-^ \faPcTPc\fo) (E3a) 

= ^JTx{TPc)P{T', a) < y/Tr(TEc) = y/l-Tr(TQc) ■ (E3c) 


In our case pid{t) = PcPid{t)Pc and Tr(rQc) = Pieak(f)- 


















